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Pre  face 


The  study  of  radioactive  fallout  from  nuclear  dctonatioj^s 
is  of  great  interest  in  this  time  of  debate  over  tlie  proiiosed 
MX  missile  field.  A  counterforce  attack  against  this  complex 
could  consist  of  thousands  of  large  yield  nuclear  weapons 
detonated  on  or  near  the  ground.  Such  an  attack  could  have 
more  far-reaching  consequences  than  most  military  planners 
consider.  A  highly  survivable,  mobile  missile  system  would 
have  little  value  to  the  millions  of  people  l.illed  froDi  the 
fallout  produced  by  such  a  counterforcc  attack. 

Presented  within  is  a  simple,  efficient  procedure  for 
accurately  determining  this  collateral  damage  of  fallout  for 
any  scenario  involving  many  bursts.  This  method  is  designed 
for  an  operational  planner  to  easily  "scope  the  pjoblem" 
without  utilizing  much  computer  time, 

I  am  grateful  to  Dr,  Charles  J.  Bridgman  for  his  guidance 
in  the  development  of  this  procedure,  and  to  ?i!y  lovely  v/ife, 
Michaela,  for  her  patience  and  support. 

John  P.  Grand ley,  Jr. 


(This  thesis  was  typed  by  Sharon  A.  nahricl) 
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Abstract 

A  method  is  developed  for  calculating  fallout  deposition 
downwind  from  a  massive  nuclear  attack  on  a  small  target 
area  over  a  short  time  span.  This  is  accomplished  using 
existing  smear  codes  and  replacing  their  existing  horizontal 
activity  distribution  with  an  approximating  function.  This 
function  is  the  difference  between  two  cumulative  normal 
functions  which  are  shown  to  result  from  superposition  of 
individual  bursts.  A  comparison  is  made  between  the  contours 
predicted  by  this  new  code  and  contours  predicted  by  the  old, 
time-consuming,  iterative  procedures.  The  new  code  has  been 
employed  in  several  different  scenarios  involving  the  proposed 
MX  field  to  determine  the  resulting  dose  contours  from  a 
massive  attack  against  that  field. 
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A  MULTIBURST  FALLOUT  MODEL 


FOR  OPERATIONAL  TYPE  STUDIES 

I .  Introduction 

A  new  Intercontinental  Ballistic  Missile  (ICBM)  system 
will  be  deployed  in  the  Western  United  States  within  the 
next  several  years.  This  new  system,  designated  the  MX, 
will  greatly  enhance  the  survivability  of  the  land-based 
leg  of  the  Triad.  This  survivability  will  be  accomplished 
through  mobility;  a  single  missile  will  be  shuttled  ran¬ 
domly  among  23  shelters  on  a  racetrack.  The  current  propo¬ 
sal  calls  for  200  missiles  to  be  purchased,  with  200  race¬ 
tracks  to  be  built  in  Utah. and  Nevada  (Ref  1:3). 

This  method  of  "hiding"  the  MX  missiles  creates  tremen¬ 
dous  problems  in  targeting  for  any  enemy.  One  possible 
targeting  option  would  be  to  attempt  to  destroy  as  many  MX 
shelters  in  as  short  a  time  as  possible.  This  option  would 
have  a  high  probability  of  destroying  a  large  percentage  of 
the  MX  missile  force,  while  interfering  with  the  command, 
control,  and  communication  necessary  to  effect  a  retaliatory 
strike.  Also  of  great  concern  is  the  corresponding  downwind 
fallout  effects  from  such  a  massive  strike.  There  are  several 
single  burst  nuclear  fallout  models  available  which  could 
predict  this  collateral  effect  by  superposition  of  individual 
burst  fallout  patterns. 


Presently  there  are  many  different  individual  burst 
models  being  used  by  many  different  government  agencies.  Mo 
of  these  models  have  their  analysis  based  on  a  study  done  by 
the  Weapons  Systems  Evaluation  Group  (WSEG-IO)  in  1959. 
However,  Bridgman  and  Bigelow  (Ref  2)  have  recently  shown 
that  the  heart  of  the  'WSEG  calculations  is  in  error.  They 
proposed  an  alternative  method  of  calculation,  which  was 
used  by  Colarco  (Ref  3)  to  construct  an  improved  single¬ 
burst  fallout  code  for  operational  type  studies.  Colarco 's 
code  produces  results  which  are  closer  approximations  than 
WSEG-10  to  the  Defense  Land  Fallout  Prediction  System 
(DELFIC) .  DELFIC  is  considered  by  many  to  be  the  benchmark 
of  fallout  computer  codes.. 

One  of  the  unfortunate  aspects  of  these  single-burst 
codes  is  that  they  cannot  easily  compute  the  effects  of 
many  bursts  within  a  relatively  small  target  area  (such 
as  an  MX  field).  Fallout  analysis  of  a  counterforce  attack 
on  a  missile  field  (or  similarly  distributed  target)  is  now 
done  by  superpositioning  hundreds  of  individual  bursts  by 
an  iterative  procedure.  Because  of  the  computer  time  needed 
for  these  iterations,  this  is  a  slow  and  expensive  process. 
Also,  because  of  computer  limitations,  this  i.iay  not  be  a 
completely  accurate  procedure. 

In  the  following  chapters,  a  method  will  be  developed 
to  produce  a  fast-running  computer  code  to  accomplish  multi¬ 
burst  calcualtions  without  the  need  for  superpositioning. 


This  code  will  be  designed  for  operational  use;  that  is, 
it  will  be  a  fast,  inexpensive  tool  for  operational  analysts 
to  use  in  predicting  a  reasonably  accurate  fallout  deposition. 
This  code  will  then  be  applied  to  a  counterforce  scenario 
against  the  MX  field  in  Utah/Nevada,  The  results  from  this 
study  will  then  be  compared  to  the  results  obtained  from 
using  WSEG-10  calculations  in  the  same  scenario. 


1 1 .  Calculation  of  the  Multiburst  lUstribution 

lVSEG-10  predicts  fallout  dose  rates  on  the  ground  with 
the  following  equation  (Ref  4:17): 

00 

6j^(x,y)  =  k  /  f(x,y,t')  g(t')  dt"  (1) 

o 

where  k  is  a  source  normalization  constant,  g(t")  is 
the  activity  deposition  rate,  and  f(x,y,t")  is  the 
normalized  horizontal  activity  function.  This  activity 
function  is  a  bivariate  normal  function  with  a  time  varying 
standard  deviation  in  the  cross  wind  direction; 

f(x,y,t)  =  -  e  •  -  6 

/2it  0  >^2tt  o  ft) 

A  / 

(2) 

where  x  is  the  downwind  distance,  y  is  the  crosswind 
distance,  t  is  time  and  a  is  the  standard  deviation. 
is  assumed  to  be  a  constant  wind  which  rotates  uniformly 
with  altitude.  This  results  in  a  constant  value  of  wind 
shear  for  the  determination  of  o^(t) 

It  can  be  seen  that  Eq  (2)  can  be  written  as 

fCx,y,t)  =  f(x,t)  •  f(y,t)  (3) 
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V.ili.  t  h  i  suhst  itutioii  it  can  I'c  si.uv.i;  ( Re  1'  a)  iRai  i  ll 
reduces  to: 


*  s(i.) 

I)(x,y)  =  k  f(y,t) 

X 


(4) 


where  t  is  arrival  lime  of  the  cloud  and 
a 


f  fy.t) 


/2r  o  (t) 

y 


(5) 


It  is  obvious  from  F.q  (S)  that  f(>’,t)  is  a  normal 
function  descriliing  the  crosswind  spread  of  a  single  nuclear 
cloud.  If  two  or  more  clouds  arc  in  close  proxijiiity  and 
merge  at  some  point  to  become  one  large  cloud,  a  single 
normal  function  will  no  longer  account  for  h.orizontal 
activity;  there  will  be  some  cumulative  effect  from  each 
contributing  single-burst  cloud.  The  fallout  on  the  ground 
from  this  new,  large  cloud  will  still  be  described  by 
liq  (4),  except  that  f(y,t)  will  no  longer  be  given  by 
Eq  (5). 

f(y,t]  for  Multiple  Bursts 

If  an  observer  were  standing  some  distance  downwind 
from  a  target  complex  (for  example,  a  missile  field!  and 
was  able  to  observe  a  large-scale  nuclear  attack  cn  t’lat 
complex,  he  would  initially  sec  many  single  nuclear  clouds 
rising  from  the  surface.  If  this  particular  attack  were 
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'  >';u'  dc  t  c'l'.a  t  i  oi'.s  ,  all  .’.oiii.;  off  v,it.h:r.  a 

sb.ort  time  of  one  anotf.er,  and  d  i  st  ri  luit  cd  nuiformly  alon” 
an  SO  mile  line  pei'pend  i  cula  i’  to  the  wind  direction,  tlien 
the  observer  would  see  100  individual  overlapping  clouds. 

If  the  liojizontal  activity  distribution  witliin  each  cloud 
alont;  tins  line  is  given  l>y  liq  (5),  tltcn  the  observer  '.vould 
sec  100  normal  functions,  as  in  Figure  1.  Note  that  y  is 
now  the  crossuind  distance  as  measured  from  the  ceiiter  of 
the  target  field. 


It  can  be  seen  from  Figure  1  tliat  any  point  downwind 
from  this  line  of  clouds  will  receive  activity  froi,)  not 
only  tlie  .single  cloud  directly  upwind  from  this  point, 
but  also  from  adjoining  clouds.  Tliis  additional  activity 
becomes  even  more  pronounced  as  the  clouds  begin  drifting 
with  the  wind  as  the  horizontal  activity  standard  deviation, 
Oy  ,  increases  with  time  (as  per  the  IvSFG-lO  analysis). 

The  total  activity  at  any  point  downwind  can  be  found  by 
finding  tlie  contr  ibut  io7i  of  each  single  cloud  and  adding 
all  100  contributions  together.  If  tlie  standard  deviation, 
Oy  ,  is  large  witli  respect  to  the  intercloud  distaiice, 
the  addition  operation  can  be  replaced  by  an  integral: 


fCy,t) 


w/2  j 

/  - e  > 

-w/2  /2',i  Oy(t) 


(6) 
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where  w/2  is  tlic  half-field  \,'idth  (in  the  example,  40 
miles)  and  N  is  the  total  number  of  bursts  (licrc,  100). 
This  function  can  be  easily  evaluated  if  it  is  split  into 
two  part  s  : 


f(y,t) 


N 

w 


-hi 


-'Ztt  a^,(t) 


>’->'0  p 


-w/2 


y-y.  ^ 

'  o 


v  2tt  o  (t ) 

y 


(7) 


which  is  easily  transformed  into 


f(y.t)  = 


•  y+w/2 

J  ^  -i. 


“-2Z  ‘ 


/  2tt 


dz 


y-w/2 

CTrr 

J  — —  e  ‘  d  z 

-00  /  211 


(8) 


by  letting 


y-yo 


(9) 


and 


dz 


-dy 
^  0 

OTT 


(10) 
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Thus,  f(y,t)  for  many  bu.rsts  is  tlic  differc-noe  ! ’(■  i  ■.■.ocn 
the  cumulativ^e  normal  functions  for  two  different  ary.umciil  s  . 
When  this  function  is  used  in  I’q  (1)  »  fl’C  total  dose  rate 
for  any  point  downv.-ind  can  be  found.  A  graph  of  this 
function  is  presented  in  I'igurc  2,  with  the  lOd  single 
cloud  distributions  added  for  comparison. 

Computer  solutions  for  dose  rates  from  f.q  (4)  would  be 
quick  and  simple  except  for  the  integrals  in  the  multi - 
burst  f(y,t)  .  Therefore,  the  following  a])prox  imat  ion 
was  used  for  the  integrals  in  I'q  (8)  (Ref  9:932): 

P(z)  =  1  •  7  (.196834-  +  .115194z2 

+  .0003442^  H  .0195272^)''  (11) 

where  z  is  the  upper  limit  of  each  integral.  The  difference 
between  the  two  integrals  was  then  multiplied  by  tlic  number 
of  detonations  and  divided  by  the  field  width  to  produce  n 
numerical  answer  for  f(y,t)  .  This  number  was  tlicn  used 
in  Eq  (4)  to  produce  dose  rates. 

As  will  be  seen  in  the  next  chapter,  dose  rates  from 
Eq  (4)  using  this  new  f(y,t)  agree  very  well  with  dose  rates 
derived  by  addition  or  superposition.  iiowever,  several  cau¬ 
tions  need  to  be  mentioned.  This  f  (y)  function  is  meant  to 
be  used  in  situations  where  there  arc  many  bursts  in  a  rela¬ 
tively  small  physical  area.  Also,  the  individual  bursts 
positioned  on  tlie  field  width  line  cannot  be  more  than  one 
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standard  deviation  apart;  that  is,  the  number  of  boi':bs 
detonated  in  a  defined  area  (the  bomb  density)  must  be 


large  enough  to 
approximate  the 
must  be  treated 


allov;  the  multiburst  f(y)  function  to 
cumulative  effect.  Smaller  bomb  densiti 
as  an  aggregation  of  single  bursts. 
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III.  Validation  of  the  Multiburst  Distribution 

Any  model  must,  by  necessity,  be  shown  to  actually 
conform  to  the  reality  it  purports  to  represent.  There  are 
several  ways  in  which  a  model  can  be  validated,  with  most 
validation  methods  comparing  model  results  to  actual  data. 

As  there  is  little  data  available  on  downwind  fallout 
dosages  from  multiple  bursts,  a  different  approach  was 
necessary.  This  involved  a  comparison  between  the  results 
of  the  new  model  and  results  derived  from  an  iterative 
superpositioning  method  of  the  same  scenario. 

Scenario 

A  scenario  involving  a  counterforce  attack  against  the 
missile  field  at  Whiteman  Air  Force  Base,  Missouri  (Ref  7:32) 
was  used  as  a  basis  for  comparison.  This  particular  missile 
field  encloses  150  silos  in  an  area  roughly  90  miles  on  a 
side.  Two  one-megaton  devices  were  then  simultaneously 
detonated  at  each  silo,  for  a  total  of  300  bursts.  These 
300  bursts  were  then  distributed  evenly  on  a  98  mile  line 
perpendicular  to  the  average  west  wind  of  20  miles  per  hour, 
and  a  shear  of  one  per  hour.  Fifty  percent  of  the  yield  of 
each  burst  was  from  fission.  The  coordinates  of  the  contour 
lines  of  the  following  unit  time  reference  dose  rates  were 
then  computed:  1000  roentgen/hour  (r/hr),  500  r/hr,  and 
100  r/hr. 
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Description  of  Superpositioning  Procedure 

A  single  burst  code  similar  to  Colarco’s  model  (Ref  3) 
was  used  for  the  iterative  procedure.  A  grid  was  established 
with  2000  miles  on  the  x  (downwind)  axis  and  182  miles  on 
the  y  (crosswind)  axis.  The  x  axis  was  broken  down  into 
201  increments,  while  the  y  axis  was  divided  into  27. 

For  more  efficient  computer  operation,  the  150  silos  were 
broken  down  into  15  groups  of  10  silos  per  group.  These 
15  groups  were  then  placed  on  grid  line  x  =  6  and  on  every 
grid  line  between  y  =  7  and  y  =  21  ,  resulting  in  15 

silo  groups  evenly  distributed  on  a  north-south  line  98 
miles  long.  To  simulate  20  bursts  at  each  silo  group,  the 
computed  single-burst  dose  rate  at  each  silo  group  was 
multiplied  by  20.  The  superpositioning  method,  then,  was 
accomplished  by  stepping  from  point-to-point  on  the  grid 
and  calculating  and  adding  the  contribution  of  dose  rate 
from  each  silo  group  to  that  point. 

Machine  limitations  dictated  the  coarseness  of  the  grid. 
Consequently,  there  was  a  reduction  in  accuracy.  A  summation 
of  all  grid  points  did  not  total  up  to  the  entire  activity 
produced  by  the  detonations;  approximately  five  percent  of 
the  activity  was  missing.  This  missing  activity  was  not 
deemed  important  to  the  comparisons,  and  is  believed  to  be 
the  result  of  the  coarse  grid  at  the  extreme  downwind  limits 
of  the  fallout  field. 
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Description  of  Multiburst  Procoduio 


Results  using  the  new  model  were  obtained  using  the 
multiburst  code  shown  in  Appendix  A.  Tlie  required  inputs 
were  inserted,  and  tlie  coordinates  for  l!ie  desired  contours 
produced.  No  machine  limitations  were  encountered  using 
this  model. 

Comparison  of  Nhil  t  ilnirst  to  Superposition 

The  comparison  v.'as  made  using  the  same  inputs:  300 
one-megaton  bursts,  each  with  50  percent  of  its  yield 
from  fission;  wind  of  270/20  willi  a  shear  of  one  iir  ; 
field  width  of  9S  miles.  Tlie  contours  obtained  were  for 
doses  of  100,  500,  and  1000  r/hr. 

Two  different  sets  of  contours  were  dra'''n  for  two 
different  particle  size-activity  distribution  curves.  The 
WSEG  contours  use  a  size-activity  curve  witli  parameters  of 
a  =  44  and  K,n3  =  .690  .  A  second  set  of  contours  was 

generated  using  s izc - act ivity  curve  parameters  of  a  -  37 
and  £.ng  =  1.5  28  .  These  parameteis  more  closely  resemble 

those  used  by  the  DliLFIC  model  (sec  Tigure  31.  .'Ml  calcu¬ 
lations  using  these  parameters  and  the  Bri  dgman/ 1',  i  re  1  ow 
procedure  for  activity  deposition  rate  (Ref  2:29-3u)  .vill 
be  termed  AFIT  calculations. 

Figures  4a  through  4c  show  tlie  three  desired  contour 
comparisons  using  the  W.SF.G  formulation  for  activity  depo¬ 
sition  rate,  while  Figures  5a  througli  5c  show  tiie  same 
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Figure  4a.  WSEG  Comparison 
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Figure  4c.  WSJiC  Comparison 
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DOWNWIND  DISTDNCE  (miu-;s) 
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comparison  using  tlie  Al'IT  activity  deposition  formulation. 

As  can  be  seen,  these  comparisons  of  superpositioning  and 
the  multiburst  code  are  in  good  agreement.  For  smaller 

dose  rates  there  is  a  difference  in  the  crosswind  deposition 
of  activity  as  dow'nwind  distance  increases  from  ground  zero. 
This  discrepancy  is  due  to  the  way  the  300  bursts  arc  distri¬ 
buted  in  the  superpositioning  procedure;  that  is,  the  15 
groups  of  20  detonations  more  closely  resemble  fifteen 
20-megaton  detonations.  This  aggregate  cloud  will  taper  off 
much  more  than  an  aggregate  cloud  of  300  one-megaton  bursts, 
as  shown  by  Figure  5c.  Consequently,  the  multiburst  contours 
are  considered  to  be  more  accurate. 
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I V  .  I-mp I  oyiiioiit  of  the  Mill  tiinn'st  Co < i e 


To  demonstrate  the  case  with  wliich  tliis  new  model  can 
be  used,  the  fallout  patterns  from  a  counterforcc  attack  cn 
an  MX  field  were  generated.  Several  different  scenarios 
were  devised,  and  tlie  appropriate  inputs  made  to  the  code 
MULTI,  which  is  listed  in  Appendix  A.  MULTI  incoiporatcs 
the  multiburst  f(y,t)  distribution  developed  in  Cl'iapter 
II.  .Xn  interesting  feature  of  this  cede  is  the  ability  to 
use  either  the  V.'SL’G  g(t)  (activity  deposition  rate)  or 
the  AFIT  g(t)  in  dose/dose  rate  computations.  Tlicrefcre, 
for  all  scenarios,  a  comparison  study  was  made  of  contours 
generated  using  the  two  different  g(t)'s 

Scenarios 

The  MX  field,  consisting  of  4600  shelters,  is  proposed 
to  be  situated  in  tlie  states  of  Utah  and  Nevada.  Several 
different  scenarios  were  created  around  these  shelters, 
based  on  several  parameters.  These  parameters  include  the 
number  of  attacking  reentry  vehicles  (RV) ,  the  yield  and 
fission  fraction  of  each  RV,  and  the  average  continental 
winds.  Two  different  dose  contours  \v-erc  generated  for  eacli 
dose  using  the  two  different  g(t)'s  .  Based  on  all  these 
variables,  a  total  of  64  plots  (or  32  comparisons)  were 
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created . 


To  simulate  a  full,  counterforce  attack,  4600  RV's 
were  targeted  against  the  MX  field.  The  crosswind  width  of 
this  field  was  estimated  at  190  miles.  The  yield  of  each  RV 
was  allowed  to  be  one  of  two  values;  either  one  megaton  or 
500  kilotons.  These  yields  are  commensurate  with  the  war¬ 
head  yields  of  Soviet  ICBM's.  A  50  percent  yield  due  to 
fission  was  assumed  for  each  warhead.  As  another  comparison, 
the  number  of  RV's  was  decreased  to  2300  with  warhead  yields 
remaining  the  same. 

Several  average  winds  were  obtained  (Ref  8)  for  the 
continental  United  States  as  a  whole  at  an  altitude  of 
40,000  feet.  All  winds  were  from  the  west  (270°)  and  had 
different  velocities  based  on  the  season.  An  average  summer 
wind  had  a  velocity  of  35  miles  per  hour  (or  written  as 
270/35),  while  an  average  winter  wind  was  given  as  270/77. 
Strong  seasonal  winds  were  derived  by  adding  one  standard 
deviation  to  the  average  winds.  This  resulted  in  a  strong 
summer  wind  of  270/70  and  a  strong  winter  wind  of  270/119. 
These  winds  were  then  applied  to  each  scenario. 

Two  different  dose  contours  w'ere  generated  for  each 
scenario,  the  doses  being  1500  and  500  rems.  Since  most 
single  story  residences  above  ground  have  a  protection  factor 
of  tKree  (Ref  7:33),  the  1500  rem  contour  represents  500  rems 
indoors.  Five  hundred  rems  is  considered  to  be  the  dosage 
necessary  to  produce  50  percent  fatalities  in  30  days 
(Ref  7:32). 
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Fatalities  will  not  be  estimated  in  the  following 
comparisons.  However,  the  contour  lines  do  enclose  the 
areas  of  high  fatalities  from  fallout.  Major  cities  down¬ 
wind  from  the  bursts  are  delineated  on  the  figures  for 
easier  reference.  It  will  be  obvious  that  the  WSEG  formu¬ 
lation  significantly  underestimates  these  areas  of  high 
fatalities . 

Results 

Sixty-four  computer  runs  were  made  interactively 
using  a  CDC  6600  computer.  Average  compilation  time  was 
1.8  seconds  with  average  execution  time  of  1.5  seconds. 

The  64  runs  were  then  combined  into  32  comparative  graphs. 

These  32  figures  were  arranged  in  four  groups  according 
to  wind  velocity.  Figures  6a  through  6h  have  a  common  wind 
of  270/35  (average  summer).  The  second  group,  Figures  7a 
through  7h,  has  a  wind  of  270/77  (average  winter),  while 
Figures  8a  through  8h  have  a  slightly  lower  wind  of  270/70 
(strong  summer).  The  strong  winter  wind  group.  Figures  9a 
through  9h,  is  last  with  a  wind  of  270/119.  All  winds  have 
a  shear  of  1/hr. 

It  is  obvious  from  all  comparisons  that  the  WSEG 
contour  lines  encompass  much  less  area  than  the  AFIT  contour 
lines.  This  is  a  direct  result  of  the  Bridgman/Bigelow 
method  of  computing  g(t)  and  of  the  parameters  used  for 
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the  particle  size-activity  distribution  curve.  Colarco 
showed  the  same  kind  of  result,  only  for  a  single-burst 
model . 


It  is  interesting  to  note  that  the  WSEG  contours  are 
probably  close  approximations  to  the  contours  generated  by 
the  Department  of  Defense  (DoD)  in  their  predictions  of 
fatalities  from  an  MX  field  attack;  the  same  basic  WSEG 
formulation  is  used.  Consequently,  the  estimation  of 
resulting  fatalities  is  not  as  high  as  it  would  be  if  the 
AFIT  contours  were  used,  especially  in  the  low  wind  situa¬ 
tions.  As  the  wind  increases  and  the  WSEG  contours  finally 
stretch  into  the  Atlantic  Ocean,  the  number  of  fatalities 
predicted  by  WSEG  and  AFIT  will  be  about  the  same  as  the 
enclosed  population  is  equal. 
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Figure  6c.  Average  Summer  Uind  Contc'urs 
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Figure  7a,  Average  U'intcr  iViiid  Contours 
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Figure  7b.  Average  V.’intcr  bind  Contour.- 
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Figure  7c.  Average  Winter  Wind  Contnu, s 
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Figure  8a.  Strong  Summer  Kind  Contours 


4.S 


LONGITUDE 


500  REM  CONTOUR.  NIND-270/70.  4600  .  Sli 


Figure  8c. 


Strong  Summer  Wind  I'.cntor,  t 


LRTITUDE  LRIITUDE 

35®  N  45* N  ^25* N  35*  N  45*  N 


RFIT  DOSE 


WSEG  DOSE 


500  REM  CONTOUR.  mND-270/70.  2300  .5MT  BURSTS. 


Figure  8d.  Strong  Summer  Kind  Contours 
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Figure  8g.  Strong  Simmer  Wind  Contour.' 
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Figure  9a,  Strong  Winter  Wind  Contours 
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Figure  Of. 


Strong  Winter  Wind  Contours 


125  W 


115  W 


105  H 


65' H 


95"  W  05"  W  : 

LONGITUDE 

1500  REM  CONTOUR.  WIND-270/119.  -iGog  .Eh:  BURSTS. 


Fij;urc  9g. 


Strong  V.’intcr  Wind  (’ontours 


LONGITUDE 


1500  REM  CONTOUR.  NIND-270/1 19.  2300  .5h 


Fif’urc  9h.  Strong  Winter  Wind  Contours 


V.  Conclusions  '  ecu:  irion t  i  oi'iS 

This  paper  has  presoT.t  eCi  a  ir.et  iru!  fcr  casii\  prcJictiiip, 
tlic  downwind  fallout  '  .r.;,  a  :'.u>.  leai'  attac!,  invol\'in'; 

man)’  dot  onat  i  oiis  on  a  t  a  !  ■  •  .  ■  v  •_  •.  i  cr  it -aria  ■■.oii.  j,  iveri 

whicii  must  he  mot  ‘  .  •  .  ■  u  ■  '  ;  1 

tact  ica  1  niu  1  t  i  '  a :  t  •  t  ■  •a;:- 

tcr i a  .  Al on ;  u  t  •  ■  .  a  •  , ;  i  t  . 

d  i  s  t  r  i  'u  u  1  .  1  I : :  ■ 
present  o  d  . 

Thi;  ‘  .n  ,  •  ,  , 

0)10  rat  i  oi  ,  ! 

The  c 0-1  o  o  ■  ' 

t i ons  oft'  ,  ■  ' 

make  f  a i t  I 
and  u i t !  .  ;  ' 

of  t  i  ir.a  and..  ■  ,  . 

were  ent  e  i'od;  .  •  ;  '  ;  '  . 

V'a  r  j  ou  :•  1  :  •  ■,  ,  , 

and  t  hi.  i  r  e'  t A  a  t  ..u  .  .  .  '  :  .  .  a  : 

des  i  r  cd  d  os  e/d  o  s  o  I  .  '  '  .  ■  .  .  ; 

tion;  wind  voloait)-  aijd  :'.i  ii  ,  mu  ■,  i  a.  :  .  ;i  !  i 

width,  i'he  line  of  Inirst.-.  i.  a!v.,m.  ■  i,  oa  to  ’  ■  •  •  ;  i  tni  ;  - 

cular  to  the  wind.  Paimuaeteis  for  tie  ;  r  t  i  a  1  i  aati\it\ 

curve  must  be  entered  if  an  A!  ’’1  calculat  I'  l,  jo  i  :  ■  .!  ;  the 

WSIiG  parameters  are  built  into  the  cc-do. 

5<d 


Recommendat ions 


The  presented  code  will  give  realistic  approximations 
of  the  downwind  deposition  of  fallout  with  great  speed  and 
ease.  This  is  exactly  what  operational  planners  need  to 
make  timely  decisions.  However,  several  enhancements  could 
be  made  to  the  program  to  increase  its  accuracy. 

The  first  enhancement  involves  wind.  The  use  of  a 
single,  constant  wind  from  the  surface  to  40,000  feet  and 
across  the  countryside  is  a  notoriously  bad  approximation. 
Some  way  must  be  devised  to  incorporate  several  different 
winds  into  the  model. 
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Ap]K‘iul  i  X  A 


l-ortr:'i:  Tfi ’>  ! 

Tliis  apjicnJi  X  contains  tlic  i  (•  1  I  ov;  i  n/,  itei.is: 

--  a  copy  of  the  code  MUI.TI  v.liich  coniputL’S 
fallout  contour;-.,  mu’,  a  ylossar}-  of  the 
teriiis  used  in  MllJ/1'1; 

--  A  user’s  guide  to  running  tiie  compiled  version 
of  tl'.c  MULTI  prograi.i; 

--  a  sample  of  output  from  M'drTl  and  .an  explan¬ 
ation  of  hov.-  to  interpret  the  output. 


MULTI  Code 

MULTI  is  a  derivative  of  SM1:.\1’,  a  single- burst  fallout 
code  written  by  Dr.  Charles  J.  Bridgman  of  the  Air  force 
Institute  of  Technology.  Several  changes  were  incorporated 
into  SMEAR  to  produce  MULTI,  most  notably  the  multiburst 
f(y)  distribution  developed  in  Cliapter  II.  MULTI  was  also 
designed  to  be  run  interactively,  with  inputs  coming  froii 
both  the  user  and  an  attached  data  tape  (TATF.d).  The  user 
inputs  are  the  parameters  of  the  particular  I'lroblem  to  be 
studied,  while  the  data  tape  consists  of  a  table  of  coef¬ 
ficients  supplied  by  Colarco  (Ref  3:67-71).  The  table  of 
coefficients  is  required  in  the  AFIT  calculation  of  g(tl  . 

The  operation  of  MULTI  is  straiglit  forward ,  with  all 
dose/dose  rate  quantities  being  calculated  as  described  in 
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Chapter  11.  The  contour  coordinates  are  given  in  terms  of 
X  (downv.'ind  distaiice  from  line  of  hursts)  aiul  y  Ceross- 
wind  distance  from  center  of  field).  These  contour  coordi¬ 
nates  are  displayed  with  the  output  at  the  terminal,  and  are 
also  written  on  a  data  tape  (TAPF/)).  T!iis  data  tape  can 
then  be  used  as  the  input  for  contour  drawing  routijies. 

The  generation  of  contour  coordinates  is  an  iterative 
procedure.  The  dow'nwind  distance  x  is  increased  by  a  small 
amount,  and  the  dose  (or  dose  rate)  at  that  x  coordinate 
and  the  center  of  the  field  (or  y  -  0)  is  then  computed 
and  compared  to  the  desired  dose.  If  tlie  computed  dose  is 
equal  to  or  less  than  this  desired  dose,  we  are  either  on 
or  outside  cf  the  desired  dose  contour  and  tlie  y  coordinate 
is  assigned  a  value  of  zero.  If  the  computed  dose  is  greater 
than  the  desired  dose,  we  are  witliin  the  desired  contour  and 
the  y  coordinate  is  increased  by  a  small  amount.  A  new 
dose  is  computed  for  this  new  (x,y)  position.  If  this 
computed  dose  is  still  greater  than  tlie  desired  dose,  tiie 
y  coordinate  is  again  increased.  This  procedure  continues 
until  the  computed  dose  equals  the  desired  dose;  1  lie  y 
coordinate  is  then  the  total  distance  incremental!;,  t  -iveled 
from  the  center  of  the  field  to  this  equality  poiiU.  Tlie 
downwind  distance  x  is  increased  again,  and  a  lu  w  y 
coordinate  is  determined  using  the  same  procedure.  Twenty- 
six  pairs  of  (x,y)  coordinates  are  generated  to  define  a 
desired  dose  contour. 
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A  glo.ssar)-  of  term?  used  uithiji  the  code  is  j/rovided 
in  Talile  A-1.  These  terms  arc  presented  in  order  of 
appearance  witliin  tlie  prograiii. 

User's  Guide  to  MULTI 

Table  A-2  sliows  how  tlic  parameters  of  a  problem  are 
entered,  Tlie  first  directive  sets  up  the  internal  switches 
for  the  program: 

M  determines  the  type  of  g(t)  calculation 
(0  =  ]VSI:G;  1  =  AT  IT); 

--  MD  determines  the  type  of  output  desired 
(0  =  dose  rate;  1  =  dose); 

N  determines  the  number  of  contours  to  be 
generated. 

The  parameters  of  the  scenario  arc  entered  next,  with 
the  desired  contour  values  entered  first.  The  yield  of  the 
attaching  weapons  (in  megatons)  and  their  fission  fraction 
are  then  requested,  as  are  the  average  wind  velocity  (in 
miles  per  liour)  and  crosswind  shear.  The  number  of  bursts 
must  be  entered  next,  along  with  the  crosswind  v.idth  of 
the  field  (in  miles). 

If  the  desired  contours  are  to  be  comj)uted  i  i:  ter],;?  uf 
dose  (that  is,  MD  =  1),  the  total  number  of  liours  wl.ich 

the  dose  is  to  be  integrated  is  entered  next.  'i  li  i  integra¬ 
tion  over  time  is  the  Way-Wigne-r  approximation  ,!'c!  S :  391 )  . 

If  a  time  of  infinity  is  desired,  zero  is  tlie  rc  lai  i  n  vl  input. 
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Finally,  if  all  calculations  are  to  be  made  using  tlie 
AFIT  g  (t )  (that  is,  M  =  1),  the  parameters  of  the  activity- 
size  distribution  curve  must  he  entered.  To  approximate  the 
DFLFIC  default,  the  parameters  given  in  Chapter  III  sliould 
be  used. 

Interpretation  of  Output  from  MUl.Tl 

Table  A-3  is  the  output  of  the  scenario  entered  in 
Table  A- 2.  A  short  summary  of  some  of  the  input  data  is 
first  given.  "AFIT  CALCULATIO.N’”  means  tlie  program  will 
utilize  the  AFIT  g(t)  in  its  computations.  The  yield, 
fission  fraction,  wind  and  wind  shear  are  printed  out,  as 
is  the  desired  dose  contour.  The  26  pairs  of  coordinate.; 
for  the  contour  are  then  listed.  Note  that  this  listing 
is  just  "half"  a  contour;  due  to  the  symmetry  of  the  cloud, 
the  other  "half"  of  the  contour  has  the  same  x  coordinates 
but  the  negative  of  the  y  coordinates. 


ooooo  o  oo  o  oooooo  no 


coMi’UTi;!;  com: 


Pf  C'Gr  i  VUL  T  I  (  I'm-''!,  l  '■-IIT,  '  5  ^  r,  w) 

:s  r-  „r::  T I c . 

TAPEC  IS  At-  DUT^L'l  1  A- E  -''■P  VALJES. 

Di>-‘EN?ir\'  if  (zr- ) ,  ) 

C0“MC‘<(  D(«),JC(;  ), Of"  ), As 

CC'iKCS  T  SI3C,  5  I  GH,  SH'  ,?C,  V  ,  i  ''' ,  T  0  ,  TP  ,  »<,  h  ,  m:" 

CQ'^HCM  NN,w,TD 

INTERVAL  SVir^HES  = 

H=:  FOP  ■JjEG,  V=  1  POP  fit’ll 

Mn=!;  FCF  3DSE  PiTE(^/pm,  MH-l  FOP  003  E(P) 

H  IS  THE  NtJi3EP  OF  CaNTD’J^  S  TD  3E  COMPUTED 

PPINT»  ,*‘E>JTEP  H(  i-  DP  1)  ,MD<«  CP  1),AN0  CDSTOJP'Sm 

READ  '  ,  C,  -fO, 

READ  VA  LMC;  OF  CO  P''0UES7E0  IN  P/HR  OR  P 

PRINT'  COKTCn,  7ALUES=  " 

Rtf',0*  ,  (C(I)  ,  I  =  l> 

READ  VlEL'i  (IN  ^tT  )  ,  FRAO  TI  0  N  OP  YIELD 
WIND  (IN  H'l'-ip)  ,  LCD  NTNO  SHEA-;  (IE 
PRINT'  t'^EvirE-  y:  ELP,FpA  CTICN,  vaCT 
E'.EAD  '  ,  YR,  PP  ,  VX  ,S-'P 
PRINT' ,“EU  EF  f*  CF  PUPS •'’3=  " 

READ-'  ,  NH 

PRINT’  ,*':'(T£P  N-S  WIDTH  OF  F:ELr  = 

Rpfl  O'  , '/ 

W---W/2. 

IF(Mr.Erm.)=PlNT'',”^‘T‘'?.'='  time  of 
IFTHD.rC.l.)  PEA?' 

SC  (SOURCE  \OPMA  L]  ’  AT:DN  :  ONSTAM)  YIELDS  R  OP  P/HP 
IF(M,EO.  .  I  S:  =  ?,  Em' R'PP"^H 
IF(M,t  0.1.  )  S:  =  2,  ■5-E  +  "P  *FF*  YH 

YIELD  DEPEND  mr  DCNS'  A  :rp  '  cr  CL'  =‘U''E? 

HC  IE  I';  <7L0-;ET 
SIGH  51, ',0  f.  IN  ^TTfjTE 

TC  IE  IN 

HC=  f-L'' G  (Y ^  (.a  r-p- (Y>'' +E  .  : )  ( L. 'iG(  f ; 

^  '2.-2) 

SIGH  =  .  1  HD/r  .  ?  ■ 

SIPPj  =  EY^(  .  ?  +f  -  ( ALOE  f  Y'U  ‘  i.  r '  '  ■')  ' 
c  TC  =  1 .  '  7‘  ( 1  2  .  -T/  •  i-c  PC/m  . . )  -  (1  .  -  ■  X-  f 

I  -PC' H'/M,' .)  ) 

C  SET  A.PHA  A  NO  BETA  PO^  W'=:  G  CALCULATIONS 
AL  =  'f  L  . 

BT  =  .09 

C  IF  THIS  IS  WS^G  :  ALOULA-IO'.,S<I^  TO  E . 

IF  (H.pr.?)  Gcr:  & 

PFl'd'-,  "FNTpy  Alp  PA  and  beta  =  •• 

REA3  »,  ALjBT 


AS  P'lSSION, 

HP-1 ) 

LEVEL,  AND  WIN"'  iH?  = 


DURATI0N=  •• 
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C  <I_C'Fr'T  T !) 

'  a  ‘ 

C  DEFIVE  ;:nSTA\TS  fO?  GfT) 

DO  >5  J-2, h 
C(  J-1)  --Co:  FFC'K,  J  ) 

25  CCsTINUE 

DO  25  J=",9 

C  (  J-2)  =CCi  FF{7i<,  j  ) 

26  CC'^riNJE 

pRin  ?5 


■c 


c 

c 

c 

c 

c 

c 

c 

c 


55  FOP  MI  (»•  /PIT  :*LCULATTCN  ♦) 

IP  THIS  IS  A  'fSEG  :  ALOMLAT  ION,  PFINT  THAT 
ec  IF  (M,  £0.-  I  =F  I'^T  53 
65  FORMAT  (*  «IESG  Ct  l.C'.JL ATI  CN  ') 


P-^INT  OUT  TH-  -i'Ali:cr<S 

PR!  '^T  7  .  ,YS  PP 
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FQ-  v.iV 

( 

vLElO  -  ■  ,‘".2, 

pp . 

FIE-"  I  O'!  ^ 

'•.ACT  lO'i  =  •  ,  -  ; .  2) 

PRIM 

->  •: 

» 

V  V, 

75 

FO.-,MA  1 

(  ' 

N[  SIJ  =  ‘jF'.l,- 

UI^■0  SHc-k 

=  S'":’  '  S-i 

FP.IsT 

C  , 

8C 

FO-  A1 

(• 

\ 0  r I  VI !  Y-"T^E 

<■  T  t  7  .. 

1 '  U  T  i  0  '  i  ~ 

uCGiiD'H L  ') 

POT  = 

1  • 

« 

PR'  NT 

;'2  5 

n,:T 

82 

FOR^'A  1 

( 1  Y 

» ■  r  4  2  1  *  A ^ 

K  A.E 

-  »,-r. 

1 AND  -^EI  A  =  ‘ 

^  F5.0 

) 

FIND  the:  TI'-'I  0-  “AXIMUH  G  (T)  AND  CALL  17  "TM** 

66  call  GMAY 

DO  12'-  i<  =  i,N' 

FIND  UDSLOFE  1 1.TD^S-.O  TION  cY  CALLING  ‘'UPHILL'’ 
CALL  JPHILL 

RETURN  -^ITH  VALU^  FD=  T'’ 

IF  (rB.SC.--  I.)  GO'FO  12. 

FIND  no-JNSLCF-  INTERSECTION  BY  VALUING  "OWHILL" 
CALL  DWHILL 

RETURN  «»ITH  VALU-  fdR  TF 
DIVIDH  HOT  LlNl  1‘rD  2^  INlFRVALS 
OT  =(TF-T6)  ',-  4 


) 


PRINT  n^LE  TITLE  FO-''  AFFROPPIATE  CALCULATION 
IF  HO,  EC,  )  ^^INT  ?;  ,  0(0 

g:  FOFM!,T(»  CjCRniNATPE  FO  R/HR  .INEM 

IF  (HD.NE,.)  ■=>RIN'  B5,  0(0 

95  FOPHUC-  rOORCINSTPF  fo^  ',^7.1,'“  R  IS-N'ZnE  LIYf*) 

CONPUT-:  AN-l  c;i*r  X,Y  VALUv,E  F'.-R  THIS  :Cl,T'  j' 

X  IS  OOWNWIf.:!  DTSTAHC"  Ff,f,N  BURSTS 
Y  IS  C?0S5W1VD  ''ISPA  NCE  FRCM  ‘JOTLI’.'E 
00  in  I  =  1, 

T  =  FL  CAP  (I-i  )*0T  +13 

X(I)  =  VY*T 

C  -IND  OOsE  OK  The  HOT  LINE  '^CR  EACH  STE^  X  (I) 

Y(I)  = 

Z  =  ABS(  (Y(T)-M>  /SIGY  (X(I)  )  ) 
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OO  0C1000  oo 


CUM.  FJNCTIO'J  S'?  CIVEN  ^’V 

ACRtMO^rr^  AsO'fECJK,  ?.•:?? 

F  =  1  . /(  2.'- ( 1.  +  .19s  "7-!  ,  1  1919^*- Z'*  1 

'  ft  =  l.-F 
'■r=F;.-F 
FY=Nt:/  (2-  ''!) 

3  =  50'  C  (T)  /';y‘  FV 

IF(  HO.  NE.  ’)  =  0»E1  (T) 

IF  00£E  3N  THE  ^37  LINE  ]S  GrTJTE^  THAN  THE  CESI?E3  DCO- 
GO  CFObSWIf-n  until  OESIRL'J  OOSE  is  INI  ERSEC  TE'1. 

IF  (n.GT.QCO)  GOTo 
Y(I)  =  P. 

GOTO  15 

OETERH^sjfiVr  ON  CF  CRDEGWl'D  DIS-TAE'CE  NEEDED 
TO  INTERSECT  DESIRE'''  DOSE 

13-' 

Hft)(W=(  W+3^  STFVfXd)  )  )  ‘  13  . 

_■■■  DO  li:  J=l,HAXW 

6=A+.l 

Z  =  AeE.(  (A-'-<)/STGY(X-(I))) 

F  =  l./(2,''(i,  +  ,19fA>‘ 74,11 94*- 
*  4-,  1 1  :  3  *.'♦*  7'»  I  4,C  !  SI  E?  » 

IF(A.LT.'n  F=i.-F 
FY  =  NFy  <2'  in  »- 
2  =  SC"  G  (T)  /7X* 

IF(  MO,  NE.'*.)  0=049I(T) 

IF  (G,C-T.j(<)  )  GO  TO  IvG 

y(r)=A 

GO  TO  ij' 


1G6 

COSTT'JUE 

IC'D 

PRINT  11'  ,  r,X  (I),Y(J) 

WRITE  (  G,  112)  ,  y  (I)  ,Y  (I) 

HE 

FORHAr  (IT.  ,2'-'.3,l) 

112 

FORMAT  (l'X,2-  1  •.C) 

111 

CONTINUE 

3RINT  HE  FIML  POTnii 
X(26t  =  VX‘ 

Y(2S»  =  C  . 

PRIM  ll'.  ,7^  ,''(25)  ,Y(EG) 
W=’ITE  <•:  ,  112  )  ,K(  25)  ,Y{26) 
12'J  CONTINUE 

STOP 
END 


(.K 


J 


SU"^RC'JT1‘JE  jHUY 

':(■■- r  (♦),:■«(-  d  ,  t'.c  )  ,,'-1 " 

C(:  Tc,  ;i  :-o, :  :c;‘S  SH'-  ,sc,  ■;x,t;.  ,t  ',tf,<,  :-',m':' 

COMMON  N\’,W,TD 

C  MAXIMUM  G{T)  15  FOUND  ''y  TFaCI^i&  G(') 

C  TO  irS  PEA<  IM  .1  HOUR  IMCRE'-EnTS 
DC  4  J:l,l  I 
T  =  F.DAT  (J-  l)*.! 

OH(  J)  =5  (T  I 

IF  ( J.E  0.1  )  GOTD  u 

IF  (DM ( J) . LT.HM ( J -1) )  GOTO  5 

4  CCuriMUE 

5  TM=T-.l 
RtTUF.N 
EMD 

SUDRCJTI  iE  JOHILL 

COMMON  D(4),DM(.  ■!)  ,  C(7  )  ,AL 

COMMON  TCf SIGG , ? IGM, SH-  ,SC , VX , T9 , TM , TF , H,MD 
COMMON  N\’,W,Tn 

C  this  5U^R0UTrN£  FIMDS  THE  CN  TH=-  HOTLINE 

C  CLOSEST  TC  THE  BURSTS  TH-  OTSIRED  CONTOUR 

C  DOSE  INTERSECTS. 

C 

C  TB=  -99.  IS  A  -LAG  TO  INDICATE  THAT  J=1l)  IS  LESS  THAH  FM 

TB  =  -99. 

C  TO=INITIAL  TIME 

TO  =  -(3,»SIG0)/VX 
IF  (M.NE.;»  TO  =  0. 

T=TQ 

C  SET  or  FOR  NS  EG  (FOF  T  <  7 ‘SIGMA) 

DT  =  .1*  ASSdC) 

OM(l)  =  . 

DO  2  J  =  E,  1.  ' 

C  SET  OT  -DR  AFI"  0®  RESET  IT  FOR  WSEG 

IFCr.GE.APS  (TD)  )  n-r=,  j 
T  =  r  4-  CT 

X  =  r*vx 
7=Mf3IGY(XI 

F  =  l.^  (2.‘  (1  .4.105^  4^-4^ 

*  +.  f  T)  34  4^74  *  3+.  )1?  '^2?'‘7'  *-4)  ■  ‘A  ) 

FA=1.-F 

FB=F4-F 

FY=SN7  (2  .  ■  M  )‘'^D 

DM  (  J»  =SC*  G(  T)  /'/X*F  Y 

lF(MD.NE.il  D“(  i)  =  El  (I)  ‘  C'-‘(  J) 

IF  (DM( J) .GE.O(K) )  GOIO  3 
IFd.GT.TM)  RFTURN 

2  CONTINUE 
R-TURN 

C  INTERPOLATE 

3  DLT  =  DT*(DM( J)-D(X) )/( CM(J)-rM(J-l)) 

TB=T-DLT 

RETURN 

END 


oooooo  o  o  ooo 


«<  ir*.  -  «kk  .  4  r  t-  r.  n  .-.  4a.v>4<.  *»i>' 


S!.nc-cuTI  .' 

co-T-icN  D(»),jf-(:  '),:("•),  A:. , '1 

CO-^MCN  Tr,  51  30,  i  IGW,  S-l-  ,vC,  7X,Tc,TH  ,TF  ,K,  H,  HD 
COHHCN  NN,W,rr 

THIS  ?Ue-OllTrs)c  FIS/TS  TM-  CN  THE!  HOTLINE; 

FARTHEST  F^^D'H  THf  rij:;cTS  WHERE  THE  DESIREO  DOGE  INFERS 


JL  =  TM/.;  ♦  2, 

DO  ♦  J  =  JL,F^^ 

T  =  “L0A1  (  J-1)  -^.5 

X  =  r^vx 

Z=MXSI GY (X ) 

F=1./(2.M  l.+  .lDf  F.'A'  7+.llEl9'<*Z*7 

*■  +,'  3  3  34ii*  L)  ^>4) 

FA=l .-F 
FB=Pa-F 

FY=^M/ ( 2. ► W)«  PR 
OH  (  J)  =5C*  G  (T)  /VX‘  FY 
IF(iD.SE.:)  0'^(  J)  ='^I(T)*D^(J) 

IFOCK)  .C-r  .O'"  (  J)  )  GOTO  E 

4  CCNTINUE 
IHTE‘‘’POLME 

5  DLT  =  .?'►  (D  (<> -D  M(  J)  )/ (  C«<  J-1) -DM  (  J)  ) 

TF  =  T  -DLT 

RETURN 

END 

4.4.4  4  4  4u..f  ;  4^.  ,.4.  44^  444  *4  »  44*,K-4~44%  4  4  4  <•  44  4  4»  4  ■■ 

FUNCTION  3(S) 

COMHCN  D(*),DM(E  0)  ,  D  {  ’  ) ,  AL  ,  ?T 

COMMON  TO,  S13G,5I3^,3H-  ,SC ,  VX  ,  T  P  ,TM  ,  T  F  ,.K,  M  ,  MD 

COMMON  N-!,W,TD 

THIS  SUpROUT[NE  CC-FUTES  G  (T)  FOF  AFIT  OR  WSEG, 

IF  (K.EO.:)  GOTO  4 
IF  (S.LT«,1)  GCrC  3 

AFIT  DOM.FUT'.riGNS 

R  IS  ^article  RADIJS  (tn  -"FTE-S) 

R=D  (1)  /  (  R»  ►  5)  +C(  2)  /  (  >»  '  4)  (3)  /  (S*'  3) 

■  ■  R=R+C(R) /(?»  5 ) +;  (5) /S  +  C  (E ) +C(? ) /SORT (S ) 


C  CONVERT  METERS  TO  MICRONS 
R  =  R-»l.  E 
C 

C  A  IS  ACTIVITY-SIZE  FJNCTION 
A  =  t  . 

P=  (ALC3(0  -  ilDG(lL))/"T 

A  =  £XP(-. 5» P» F)  /  (SORT { G,2E 3) ' GT* R) 
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C  0'?/DT  IS  OSSH/fiTIV'  3^-  R 

DrOT  =  -?  (1)  /  (S  )  -  4.  ^C(2)  /  (S'  '5)  -  3.^C(  3)  /  (5 

0;?DT  =  0Srr-E  .  '•C(0/(3»‘')-C(:')/(S»S)-.S*-C(7)/(;‘*l.?) 
DPDT  =L''  DT*  1  .  F 6 
G=  A'ABS(3R3T) 

Rf  TUf.N 
C 

3  G—  -  • 

RETURN 

WSEG  COhPUT AT  ION’S 

4  X-7X‘S 
n=ABS(X) 

H=1  .  +  .19-:.3  5:<  '  0  +  .  llSlqiv'  O'^Q+.U  17^4^  3  +  ,  '  4 

PHIri.  -  l./CS.”  H**^'*) 

IF  (X.LT.j.)  PHI  =  I.-cht 
G  =  PHI*  S<=  (  -S/r  O/T'' 

RETURN 
END 

«w>«4«V44>'  +  > 

FUNCTION  31(5) 

CO'-'KON  0(4)  ,0>!  (i  .'))  »  0(7  ),AL  ,PT 
CO^fCN  TO,  SI30,  SIS'-',  SNR  , SC,  VX,T?,T>!  ,TF,t<,  M,MO 
COMMON  NN,W,rD 

THIS  FUNCTIjN  CON'/IRTS  TOSE  KATE  TO  DOSE 
USING  THE  V4VY-MIGN-R  A OXI  MAT  I OU 

TA  =  S 

IF(TA,LT, , 1) T A  =  .  1 
TE  =  TA4TD 
BI=E./(Ta»»,2) 

IFdE.GT.TA)  3I=qi-S,/(T  E"-*.  2) 

RETURN 
END 

Q44l»44'44>4r4>44  4'»  »»4444*.4,4*.#V4*» 

function  5IGY(X) 

COMMON  0(4),D»'(:.  T  )  ,  S  ( '’  ) ,  A  L.  ,  E  T 

COMMON  TO,  SlSO,SIGM,c:Hr  ,sc,  VYjTC'  ,  ,  .‘I,’'''. 

COMMON  N  ^J,TO 

C  THIS  7UKCTI0'(  DOMPJTE'^  SIS^A  Y,  AS  F£E  W3lG-1. 

C 

TS=X/VX 

IFdS.ST  .3.  )  TS  =  3  , 

TR=l.+8. ’TS/TC 

SIGY  =  S'^RT  (SIGD'  *2*  TR  ♦  (SI  GH*  SHR-*  X77  X )  *  •  2  5 

RETUF  N 

END 


7] 
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FU\'CTIOh  :DS:=^FC^K,M) 

DIKENSICN  ft  (21:  !  ,  9) 

THIS  -UNCTIJ')  RFftCL-  DATA  VAPF  A),D  IK'T  ERPOL  ATES 

FOR  THE  CCi'FIOIEVTS  ^'EEDED  FOR  THE  GIVrIN  ftLTIT'JOE  ■^< 

(REWIND  t 
00  1 t  J =1,250 

RE;  D(v,l  :  )  (A  (3  ,  J)  ,  J  =  1 , 

REA nc; ,  1  1)  (ft (I, j), j=o,  o) 

ICE  '^0-HftT  (F-  .  1,  ■•.  Ell  .5) 

IDl  FC'Pf^tT  (F  ,.  1,  3EU  .9) 

IF(7t<.LE,ft  (I,  1))  GO  TO  11 

10  CC'ITINUE 

11  X2=A(I,1) 

Xi=A(I-l,l) 

Y?  =  A(r,N) 

Y1  =  A  (:-!,'() 

C0£FF=((Y2-Y1)/(X2-X1))  *7K+ ( (  X2*  Y1  -  X 1 '  Y2)  /(X2-X1)  » 

RETURN 

EKD 
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TABLl:'  A- 1 


G 1  o s ary  of  Tcri:is  V.’itliin  !!U!T  1 
INPUT  PARAMirrr.RS 


D(I) 

Contour  value  in  dose 

or  dose  rate 

YM 

Yield  of  the  weapons 

(in  megatons) 

FF 

Fraction  of  yield  due 

to  fission 

VX 

Velocity  of  the  wind 

(in  miles  per  hour) 

SIIR 

Crosswind  shear  (in  hours 

NN 

Number  of  bursts 

W 

Width  of  the  field  (in  miles) 

TD 

Time  of  duration  for 

dose  computat i ons 

(in  hours) 

AL 

Alpha  parameter  from 

act iv ity - s  ize 

distribution  curve 

BT 

Beta  parameter  from  activity-size 

distribution  curve 

WSEG-10  PAR.A'IETERS 


HC 

Nuclear  cloud  height 

(in  kilofoct) 

SIGH 

Cloud  thickness  parameter 

SIGO 

Initial  cloud  radius 

parameter 

TC 

Time  constant 

SC 

Source  normalization 

const  ant 
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TADLi:  A-l  rCoi.tM'! 


roM^'UTATiox  rAR.'A;!;!;;:;;- 


C(J) 


MAXW 


Nuclear  cloud  hei  r, lit  [in  kilo:r,cters ) 
Coefficients  for  AFIT  e,(t)  computation 
Interval  along  liot  line 
Time  of  maximum  g(t) 

First  time  of  occurrence  of  desired 


dose  on  hot  line 


Final  time  of  occurrence  of  desired 


dose  on  hot  line 


Time  (in  hours) 

Dovv'nwind  distance  (in  miles) 
Crosswind  distance  (in  miles) 
Cumulative  normal  function 
Crosswind  distribution  of  activity 


Dose  or  dose  rate 


Crosswind  mileage  counter 
Crossv\'ind  limit  of  commit  at  ion 


FUNCTION  G(S)  rA:A'.''^-TFRS 


DRDT 


L 


Particle  radius  (in  meters) 
Activity-size  function 

Change  in  particle  radius  with  resr  .’Ct 
to  time 

Cumulative  normal  function 


TABLH  A- 2 


Interactive  Inputs 

ENTER  M  (0  or  1),  MD  (0  or  1),  AND  N  (•:/  CONTOURS)  =  1,1,1 
ENTER  CONTOUR  VALUES  =  1500 

ENTER  YIELD,  FRACTION,  WIND  LEVEL,  AND  LIND  S!1R  =  ],.5,35. 
ENTER  ^  OF  BURSTS  =  4600 

ENTER  N-S  WIDTH  OF  FIELD  =  190 

ENTER  TIME  OF  DURATION  =  0 
ENTER  ALPHA  AND  BETA  =  37,1. 5?R 
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TAKl  };  A- 3 


Sample  Oil 1 1" 1 1 1 


Vita 


John  !■ .  Crandlcy,  Jr.  was  born  on  25  November  1951  in 
Pittsburgh,  Pennsylvania.  lie  graduated  from  Milford  High 
Scliool  in  Milford,  Connecticut  in  1969,  and  attended  the 
U.S.  Air  Po.rce  Academy  where  he  received  the  degree  of 
Bachelor  of  Science  (Chemistry)  in  June  1973.  After  grad¬ 
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